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ABSTRACT 

The theoretical model accounts for drift and cyclotron components of 
ion motion in a partially ionized plasma. Density and velocity distribu- 
tions are systematically prescribed. The flux into the neutral particle 
analyzer (NPA) from this plasma is determined by summing over all charge 
exchange neutrals in phase space which are directed into apertures. Cal- 
culation of the process is continued through the NPA using appropriate 
cross section data to obtain analyzer output distributions. Theoretical 
results were compared with NPA measurements on four plasma heating de- 
vices having radial electric, E, fields and axial magnetic, B, fields. 
Drift velocity, in the azimuthal direction, is identified with E/B. Es- 
pecially detailed data, obtained by sweeping the line of sight of the 
apertures across the plasma of the NASA Lewis HIP-1 burnout device, are 
presented. Selection of randomized cyclotron velocity distributions 
about mean azimuthal drift yield energy distributions which compare well 
with experiment. Use of data obtained with a bending magnet on the NPA 
showed that separation between energy distribution curves of various mass 
species correlate well with a drift divided by mean cyclotron energy 
parameter of the theory. Use of the guiding center model in conjunction 
with NPA scans across the plasma aid in estimates of ion density and E 
field variation with plasma radius. 


INTRODUCTION 

Determination of ion motion in a plasma from diagnosis of the escap- 
ing charge exchange neutrals is made difficult by the many processes in- 
volved. Not all of the particles admitted into the apertures of a neu- 
tral particle analyzer (fig. 1) influence the output signal. Only those 
which are suitably charged after passing through a stripping cell and 
which have a resulting energy, within a small spread, are able to reach 
the output detector. This band of energy is dependent on beam resolution 
and the setting of the electrostatic energy selector. Relating the mo- 
tion of such terminal subgroups of particles back to the characteristics 
of the total ensemble of ions in the plasma involves unwieldy kinetic de- 
scriptions in phase space. Simply the selection and determination of the 
many concomitant multivalued roots in such an inverse type of problem can 
be quite tedious . 

As an alternate approach, an attempt is made to model the essential 
features of the ion motion in the plasma; A guiding center model is 
hypothesized for the crossed E and B field type of heating device of 
interest herein. Such a forward method of solution also simplifies cal- 
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culation of the processes inside the NFA. Comparison of these results 
with experimentally determined NPA output can be used to evaluate, modify, 
and extend the plasma model. 

The guiding center model gives results which agree with experimental 
data from four different plasma heating devices. This model and its eval- 
uation are discussed herein. 


PLASMA MODEL 


The plasma is assumed to be axisymmetric and have a radially inward 
directed electric field, E, normal to the axis. The magnetic field (flux 
density) B, is assumed to be in the axial direction and uniform in a 
plane of diagnostics normal to the axis. In figure 2(a) the crossed 
field velocity, v c , of a guiding center drifts clockwise about the axis 
of the plasma while the ion cyclotron velocity, v C y, is directed counter- 
clockwise about its guiding center. The small change of direction which 
a heavy particle undergoes during a charge exchange event, as well as 
the small time interval of interaction for the event to occur are assumed 
negligible. Thus in the analysis, for a particle to enter the aperture 
it is necessary that the vector sum of the ion drift and cyclotron ve- 
locities be directed toward the aperture at the time of the charge ex- 
change event. From figure 2(a) 


V cy “ v l + y c + 2v a v c sin 6 (1 > 

2 

Here v a '* 2z/m where £ is the energy of the entering neutrals which 
is set by the NPA. In general, v c is dependent on the guiding center 
location, ? c , which is displaced from the ion location by the cyclotron 
radius, ? C y. From figure 2(a) and some trigonometric manipulation* 

r sin 6 = y + r hil - (v cos 0/v ) 2 
c J cy \! c cy' 

and using equation (1) this becomes 

r sin 6 = y + r (v + v sin 0)/v (2) 

c J cy a c cy 

From figure 2(a) 


r c sin 9 = y ~ r cy sin(0 + 6) , 

where 0+0 is the angle that r C y makes with the positive x axis. 
By law of sines and relation between inverse trigonometric functions 


0 + 0 = cos ^(v cos 0/v ) 

c cy 


= sin" 1 


+ 

L 


P 


(v cos 0/v )' . 

c cy ji 
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Here the plus sign corresponds to the case where the x component of 
v c is directed away from the aperture. This second solution is possible 
only in the lower half plane. The ions corresponding to this sign have a 
guiding center region nearer to the axis than those identified with the 
minus sign. 

The magnitude of the cyclotron velocity is considered to be distribr- 
uted about some mean value. For full randomization of this constituent 
there exists a cyclotron temperature, Tcy • The sizes of the domain of 
r c , shown by the shaded area on figure 2(a), are thus due to the spread 
in v C y as well as to the aperture height. 

The number of particles entering the apertures per unit time is 
equal to that of the various species which in turn can be expressed as 


N = a 

a 



o , v n.(r )n (r )f (v)r dr d9 dv 
ch a i c n c v c c 


(3) 


Integration is over all coordinate (r c ,6) and three-dimensional velocity 
space, v, from which contributions can be made. It is assumed that the 
velocity distribution can be broken into components parallel and perpen- 
dicular to the plasma axis. 


Since the apertures have finite width, £, in the axial direction, 
there is a dependence of N a on v| | . For a randomized distribution of 



f | j ( v | | ) = v'aV (2e ^ ) 

where the constant, a, can be identified as m/2kTi t if so desired. The 
contribution of vri to equation (3) is 


(/a/2) 


~y a tan a 2 

-av j | . 

e 1 1 dvi | - (/tt/2) erf (va tan a) 


'-v a tan a 


* /a v tan a 
a 


(A) 


The quantity tan a is a geometrical constant as shown in figure 2(b) 
and the final approximation is good for small a. 

At any specific location in the plasma, cyclotron velocity can vary 
whereas drift velocity is fixed. Thus 
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f jA )dv X = f Ay + ’c )d V 


= f ' (v )dv 
cy' cy' cy 


= (1/4) v sin d) f (v Jdv dd> 
v 7 ' cy Y cy cy cy r 

r 

where <\> is azimuthal angle and noting that f C y is isotropic xn v C y 
space. Assuming a fully randomized distribution of the magnitude of 


f 

cy 


-mv 

(m/4kT cy )e 


2 

cy 


/ 2kT C y 


(5) 


Substituting equation (4) and (5) into (3) , integrating over 
ing (1), and converting from velocity to energy space gives 



, ’ch n i (t c /r D ) W r D> 


({>, utiliz- 


x exp 



2V & 
c 


sin 0 


»)<£(! 


+ V c Sin 0/>^) (r c /r D )d(r c /r D )d0 d£ 


where 


£ = mv 2 / 2kT = e/kT 

a cy cy 

C = &/a r tn a kT /2m 
D cy 

Velocity and radial distance have been made dimensionless by arbitrarily 
dividing by ( 2 kT C y/m)l /2 and ion diffusion radius, r D , respectively. 

The electrostatic plates of the NPA select small ranges, Ae/2 about 
e so that at the NPA output . 


*. • 


r 9 (0) /r T 


1 = § ■ 



WWW 


' 6 1 V ^ 0 ,-l <9)/r D 


x exp - (j£ + + 2V yffi; sin o\(l + V sin 0/v^)(r /r n )d(r /r n )d6 

( 6 ) 

where is the fractional ionization of the neutrals as a result of the 
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stripping cell,.^ is a factor (proportional to e) accounting for the 
beam resolution , and k accounts for the dependence of the output de- 
tector on the energy of the impinging particles. 

Equation (6) can be numerically integrated for arbitrary density, 
drift, and temperature distributions. It is expedient, however, to first 
hypothesize various elementary reference models in order to determine the 
sensitivity of results to the large number of free variables and to make 
preliminary comparison with experimental NPA output data. 


Two drift velocity distributions are especially expedient: (1) v c 

independent of r c /rp and (2) v c proportional to r c /rp. Case 1 may 
be argued to apply, for example, when the bulk of the plasma is outside 
of sheath regions and the applied E field is dominant. If a field in- 
duced by ambipolar diffusion, however, were dominant, then 

E = (kT e /en i )dn i /dr c (7) 

(ref. 1) . Use of a typical diffusion profile 


n i = n i<°) e 


-(r c /r D )' 


in equation (7) yields a result favorable to Case 2. 


( 8 ) 


Simulation of the applied E field by that formed with a line source 
of electrons on the axis suggests that drift varies as l/r c * Preliminary 
estimates, however, with use of this assumption did not. yield good approx- 
imation to profiles obtained experimentally by verticle scans across the 
plasma . 

For nonuniform drift the cyclotron radius must be less than some 
characteristic length over which the electric field undergoes appreciable 
change in order for the orbits, relative to v c , to be approximated by 
circles . 


Uniform Drift Velocity 


For constant T cy the integration of equation (6) over r c /rp is 
now dependent only on n^ and n n „ For an exponential density profile 
of neutrals as 


n 


n 


bn(r c /ro) 2 

n„(0)e 


and a density profile for ions as in equation (8) it follows that 
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dN. n?2 

dgf - C^JU ± (0)a n (Q)da ch /2a - b n > / (1 + V c sin &//§) 


x exp -Qt + v + 2V /if sin 6) 
c c 


x ^*p£a - V«cyv3 ■ exp [ _(i ■ v (r c,2 /r D> 2 j 


►de (io) 


where from equation (2) 


and 


= I +ji(v^ + V sin 9) J / sin 0 

D l r D C J / 


Vi = r / (r~V ) = x /2mkT '/ (Br _e) 
cy v D cy' V cy v D 


( 11 ) 


( 12 ) 


To determine r Cj 2^ r D replace y in equation (11) with y^, plus aper- 
ture height, h« The parameter, p can be considered as a mean cyclotron 
radius divided by the diffusion radius. 

The limits of integration on 6 are symmetrical about the y axis. 
They can be found simply by observing where the contribution to the flux 
reduces to zero during a numerical integration of equation (1). Approxi- 
mate limits can be first obtained by use of figure 2(a) to make this 
process more efficient. 


Drift Proportional to Radial Distance to Guiding Center 


For this E field, and constant T C y, it follows that V c is equal 
to V c (l)r c /rj).* Here V c (l) is the value of V c where r c = r . Equa- 
tion (2) then becomes 


In this case centrifugal force, and azimuthal current due to unequal 
azimuthal drifts of electrons and ions (ref. 2) can be included in v Q 
since v c remains proportional to r c . Solution of the force equation 
(ref. 3) for v c 


“ mv y r c = q( ' E + v c B) 


by completing the square gives 


'c ■ r c(* V 




“c,i' 2 ] ' 


r v (1) 
c c 
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(r c sin e)/r D - (y/r D f pv®)/ [1 ± yVJl)] £ Y <13) 

Note that the boundaries of the guiding center regions are now straight 
lines parallel to the x axis . The flux equation can now be written as 



de dY 

exp (ip /sin 0) 

(14) 


x = 1 - b 

n 

and 

n = C^|^a ch n 1 (0)n n (0) 

With the substitution £ = i|//tan 0, the integration over 0 can be per- 
formed to yield 

_ _t[j^ 

-/ir ijje v [erf(ip cot 82 ) - erf (41 cot 0-^)]/2 

, — ^2 

which is approximately equal to -/it ip/e v . The resulting final integral, 
over Y, can be simplified by transforming to integration over Z where 

Z - v® V c (l)Y +5V^(1 )/[k + V 2 (l)] 


This gives the final result 


d» - exp(^ - 2V e (l)^Y 2 )J 

d£ 2^ 3/>2 l e xp(?) 


+ ‘# 





erf 

IV-V“ + ^ s 2_|- erf 

LV V + ^ Y iJ 


exp 


( 15 ) 


*2/ = K + v~(i), 


where 
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Y 1 = ^ y l^ r D * * VV C (1)] 

and for Y£ replace y-^ in equation (16) with y^ + h. 


( 16 ) 


Numerical . Calculations 

« 

The expressions for dN/d<f (eq. (6) , (10), and (15)) are dependent 
on the parameters T C y, p, V c (l), y^/t^j ant ^ h/rj) f° r e set fc - e 
NPA. Selection of the radial distribution of V c , n^, and n Q must 
also be made. All curves will be normalized such that their peak values 
reach unity. Thus, the constants n^(0) and n n (0) of equations (8) 
and (9) do not influence results in the paper. The detector character- 
istic, Of, is assumed constant and assumed proportional to e. The 
energy distribution can then be calculated provided sufficient data are 
available to determine and 5^. 

Curves were fitted to the charge exchange and fractional ionization 
data for the species of interest (figs. 3 and 4) . When required, numer- 
ical integration was performed by use of Gaussian quadratures with an 
IBM 7094 digital computer. 


COMPARISON OF THEORETICAL RESULTS WITH NPA OUTPUT 

Theoretical results were compared with NPA output data from four 
crossed-field type plasma heating devices (refs. 4, 5, 6, and 7). Avail- 
able from all of these investigations are curves of dN/de against e 
for helium and/or various hydrogen species. In addition, traverses in 
the y direction of the diagnostics plane were made in the Lewis HIP 
and SUMMA plasma heating devices. Such vertical scans are most helpful 
in the study of radial density distributions. The energy traces, how- 
ever, were quite insensitive to density distributions, but were of prime 
importance in the estimation of T cy and v . Unless otherwise stated t 
uniform neutral (b n =0) and diffusion type ion (eq. (8)) density pro- 
files were used in the numerical calculations presented for the energy 
distribution. Bending magnets were used on the NPA of references 4, 5, 
and 7 so that only particles having a specified mass-to-charge ratio 
could reach the detector. 


Energy Distribution Parameters Ae and Ep 

The various energy distribution curves can be quite well .identified 
by the combination of their half widths, Ae, and their most probable 
energy, e p ; that is, the energy at which the NPA output current is a 
maximum. Plots of Ac against e p enable an insight into the extent, 
of influence of a number of variables on the energy distributions. Re- 
sults for He + ions charge-exchanging on He neutrals, for a B field of 
2 tesla, are shown in figures 5(a) and (b) . 
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The difference of results obtained by use of the two azimuthal drift 
assumptions is very pronounced. Results based on case 1 (V c independent 
of r c ) form a grid of widely separated lines of constant V c , whereas 
results based on case 2 (V c proportional to r c ) form such a narrow band 
that they can be represented as one straight line which is independent of 
B as well. For hydrogen this band is slightly curved and not quite so 
narrow . 

Constant y lines coincide with those of constant T C y; see equa~ 
tion (12). Where y increases beyond unity, the mean cyclotron radius 
is greater than the diffusion radius and use of a guiding center model 
becomes progressively less justifiable. 

Values of Ac and £p obtained from NPA otuput traces on an x-y 
recorder are shown by the data points on figure 5 (b) . Two operating 
lines can be distinguished. One is for the cylindrical type anode used 
in both HIP and SUMMA; the other for a water-cooled ring anode used only 
in SUMMA. Both lines are to the right of the case 1 line and well into 
the region of case 2. 

Drift velocity, v c , is next identified with E/B to show lines of 
various E values for constant B in figure 5(b). Quite reasonable 
values of E and T cy were found in the interpretation of NPA results 
by use of such plots as long as y was less than unity. These variables 
are more descriptive than V c and y and will therefore be presented in 
the following discussion. 


Vertical Scans 

A typical curve of current from the NPA output detector against 
y/r D is shown by the line comprised of short dashes in figure 6(a). 
This trace from HIP was obtained with the electrostatic plates (ref, 8) 
at a fixed energy setting of 3.5 keVi This is near e p for the test 
conditions which were 13.5 kilovolts across the electrodes and a 
1.1 tesla B field. 

Vertical scan results are shown for both the constant drift and 
the drift proportional to r c theory even though this latter theory 
does not give energy distributions which agree with NPA data at high 
V c (fig. 5(a)). 

The shape of the energy distribution curve (not shown) was well 
approximated by assuming Tcy - 0 .65 keV, h/rj) = 1.6, and a uniform 
E = 2700 V/cm. The visible radius of the plasma was approximately 
2 centimeters in HIP, thus rp was set at 1 centimeter, (With a • / .• 
diffusion’' . profile • approximately : 98 percent iv of the plasma lies 
within r c /rp =2.) Using these conditions, with tha exception that 
h/rp = 0.8, in the V c « r c theory gave results which also agree quite 
well with the shape of (although vertically displaced from) the NPA 
trace. The absolute y/r D location of the experimental curve is not 
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known within about 0.5 centimeter, so that only the shape and relative 
y/r D location of the theoretical curves can be appraised. The theoretr- 
ical curves have an asymmetry about their peaks which is not as apparent 
in the experimental curve. Figure 6(b) illustrates the sensitivity of 
the profile to h/r D . Results for both E type of drift variations are 
insensitive to aperture height below an h/rp of 0.2. 

At high values of T cy the density of neutrals should be consider- 
ably reduced at the center of the plasma. An arbitrary density distribu- 
tion of neutrals using b n = 0.75 was calculated in conjunction with the 
previous conditions . This profile has relatively few neutrals at the 
plasma axis (5 percent) compared to that at the outer edges (where 
r/rj) = 2) consistent with trends of references 9 and 10. The shape of 
the profile (fig. 5(a)) was quite insensitive to changes in b n , but the 
peak of the curve shifted to higher y/r^ with increase of b n . Locar- 
tion of the profile peak may thus be useful in determining neutral density 
profile. 


Energy Distribution Curves 

A closer determination of E and T cy values from the NPA traces, 
than those obtained from Ae and e m p values, can be made by comparing 
entire curves. Figures such as 5(a), however, provide excellent starting 
points from which the E and T cy values can be adjusted for better 
curve fits. 

Figure 7(a) shows a family of curves obtained from SUMMA operating 
with helium. The B field was varied from 1,1 to 3.3 tesla while the 
voltage (across a cathode-anode combination at each end of the plasma) 
was held constant at 5 kilovolts. The theoretical lines were obtained 
using the uniform drift theory with h/rp = 1.6. Note that the whole 
family of NPA traces can be quite well matched with consistent guiding 
center model curves . 

Plots of E, V c , y, and T values, determined by the theory, 

against B is shown in figure 8. Results indicate that increasing B 
enables maintenance of higher radial E fields in the measuring plane 
midway between the two sets of electrodes . Increasing B reduces the 
drift energy. Since drift energy is in part transferred into random 
energy it is quite plausible that T cy should decrease with increase of 
B. At low B .values but at higher voltages it was found that can 

be considerably less than one. Thus the drift energy can.be much less 
than the random energy indicating contributions also from another mech- 
anism, such as resistive and/or turbulent heating. 

Applications of the guiding model to data obtained from the Lewis 
Bumpy Torus operating with deuterium is shown in figure 9. The data 
presented here, although nontypical of this device, were selected since 
B and T cy are especially low causing V c to be very high. Theory 
predicts that the separation between the e p of different mass-charge 
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ratio species increases with increase of V c .. 

The curves. of both D + and D 2 were quite well fitted with 
E = 1 kV/cm and T cv =100 eV in the theory. Relatively flat ion 

y -r /r 

density profiles (n^ = n^(0)e c ®) were used, guided by the results 
in reference 6 for the high current mode. 

Data from the Oak Ridge Burnout V apparatus are shown in figure 10 
for atomic hydrogen (ref . 7). It is quite significant that the value 
of ion temperature reported in this reference is very close to the value 
CT C y.= 1.25 keV) giving the best guiding center model agreement with the 
NPA data. An E field of 6,5 kV/cm was obtained from the theory. This 
is reasonable for the 10 kV applied at the electrodes in the experiment. 


CONCLUSIONS 

The guiding center model offers explanation to a number of the domi- 
nant features of NPA output traces obtained from crossed-field heating 
devices which have a characteristic azimuthal drift. The half widths as 
well as the most probable energy of the NPA energy distributions were 
found to increase with temperature of the cyclotron component, with mass- 
charge ratio of the ions, and with the strength of the radial electric 
field. The separation of output peaks between two different mass-charge 
species in the same plasma increases with the ratio of drift velocity to 
mean cyclotron velocity. Of the two radial drift distributions studied, 
the one independent of radial distance from the center of the plasma 
yielded results which best fit the NPA data. 

The guiding center model indicates that NPA energy distributions 
are quite insensitive to radial density distributions of both ions and 
neutral background particles. Vertical scans obtained in HIP and SUMMA 
heating devices with the NPA energy selector held constant gave contours 
of current against vertical displacement across the plasmas which agree 
well with guiding center model results based on a diffusion density pro- 
file for ions. Theory indicates a vertical shift of the curve peaks with 
decrease of neutral density near the plasma axis of symmetry. Increase 
of drift velocity causes a shift in the opposite direction. 
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APPENDIX - SYMBOLS 

[The International System of Units is used unless otherwise stated.] 
A mass number 

a constant in definition of parallel velocity distribution function 
B magnetic flux density 

b constant in equation (9) 

C constant in equation (6) 

E electric field strength 

<£ dimensionless energy, .e/kT 
£fr stripping cell efficiency 
f velocity distribution function 
beam resolution 
h aperture height 

detector efficiency 
k Boltzmann constant 
% aperture length 

m mass 

• 

N number of particles per unit time 

n number density 

q charge 

r radius 

T temperature 
V dimensionless velocity 

< + vJ(d 

v velocity 
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x,y,z Cartesian coordinate; x and y are in the diagnostics plane; 
z is parallel to the plasm axis 

Y defined by equation : (13) or (16) 

Z space used in integration of equation (14) 

-y —y 

a angle between vy and. v^. 

3 angle between cyclotron and guiding center radii 

e energy 

n coefficient in equation (14) 

0 angle between cyclotron radius and x axis 

k 1 - b 

n 

p ratio of mean cyclotron radius to diffusion radius 

5 $/tan 0 

a cross section 

<j> azimuthal angle about B in velocity space 

I l> + v*(l) Y 

ca cyclotron frequency 

Subscripts: 
a aperture 

c guiding center 

ch charge exchange 

cy cyclotron 

D diffusion 

e electron 

1 ion 

n neutral 

p most probable 
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v velocity 
|| parallel 
-1 perpendicular 
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Figure 1. - Neutral particle ana iyzer. 
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(a I NEUTRALS ENTERING APERTURES IN x,y PLANE NORMALTO PLASMA AXIS. TWO VERTICAL 
POSITIONS OF THE APERTURES ARE SHOWN. 

Figure 2. - Guiding center model vector diagrams. 
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(b) NEUTRALS ENTERING APERTURES IN x,Z PLANE PARALLEL TO 
PLASMA AXIS. 


Figure 2. - Concluded. 
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Figure 4. - Fractional ionization in stripping celt. Data of reference 21 
extrapolated to 0. 1 torr. 



(a> USE OF V c AND u GRIDLINES. EXPERIMENTAL DATA OBTAINED AT B FIELDS 
RANGING FROM 1.0 TO 3.5 TESLA. 

Figure 5. - Uniform E performance grid for helium. 6 55 2 tesla, hl'r D = 1.6, b n = ( 
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{b) USE OF E AND T cy GRIDLINES. 


Figure 5. - Concluded. 



Figure 6. - Comparison of theory with helium data from HIP. Anode 
voltage set at 13. 5 kV. B * 1.1 tesla; E at r c /r D = 1 set at 2700 V/cm 
in theory; T C y - 0. 65 keV. 



lb) EFFECT OF APERTURE HEIGHT, E « r c /r D . 
Figure 6. - Concluded. 
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ENERGY, £, keV cs-?ooa 6 

Figure 7. - Comparison of uniform E theory with measurement in 
SUMMA using helium. Electrode voltage = 5 kV, h/r D - 1.6. 
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ENERGY, e. keV cs-70084 

Figure 9. - Comparison of theory with data from Lewis 
Bumpy Torus using deuterium. Anode voltage set at 
5 kV. B * 0 . 2 tesla; h/r D - 0. 14 ; n f * expiry r D >; 

V c /VA a 3.6; |j/VA= G.65. 


3 TEMPERATURE, RADIAL ELECTRIC 
NPA OUTPUT, dN M e S T cy- keV F,ELD ' E - kv/CM 
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Figure 8. - Results determined from guiding 
center model with curve fits shown in fig- 
ure 7. 
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ENERGY, E, keV cs-7 ooei 

Figure 10. - Comparison of theory with Oak Ridge 
Burnout V data for H + charge exchanging on H. 

B » 2.5 tesla; hfr D - i ; uniform E - 5 kV/CM ; 

T cy - 1 keV. 



